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We consider a deformable body immersed in an incompressible liquid that 
is randomly stirred. Sticking to physical situations in which the body departs 
only slightly from its spherical shape, we calculate the diffusion constant of 
the body. We give explicitly the dependence of the diffusion constant on the 
velocity correlations in the liquid and on the size of the body. We emphasize 
the particular case in which the random velocity field follows from thermal 
agitation. 
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Deformable bodies that are nearly spherical can be found in various soft matter systems 
such as emulsions and complex fluids, where deformable bodies are immersed in a host fluid. 
It is easy to envisage a situations in which the host fluid is stirred randomly For 
example, mixing of bodies in an emulsion by external mechanical vibrations (e.g., ultra- 
sound). Another example is thermal agitation of a complex fluid. Hence, it is interesting to 
investigate the dynamics of deformable bodies in random flow. In this work we study the 
movement of a single deformable body in a random velocity field, which is uncorrelated in 
time and correlated in space in a general way. We focus on three aspects of this subject: 
derivation of the explicit equation of motion for the body's center, the Mean Squared Dis- 
placement (MSD) of the center and the effect of the body's size on the diffusion constant. 
Consider a single deformable body immersed in a host fluid. The system is chosen to have 



the following characteristics: 



1. The host fluid is incompressible, V • v — 0. Moreover, we assume that the Reynolds 



number is small so that the Stokes approximation is applicable. 

2. The body is characterized by an energy that depends on its shape. The shape of 
minimum energy is a sphere. The energy may be surface tension |§, Helfrich bending 
energy ||||, etc. The body's surface shape is described by the equation: ip{f) = 
where ip(r) is a scalar three dimensional field. 

3. A deformation of the body induces a force density on the host fluid. As a result, the 
fluid's velocity is given by the sum of v ex t, which is caused by external sources, and 
Vip, which is induced by the body. 

4. The external velocity, v ex t, is random and is chosen to have zero average and known 
correlations. It is convenient to define the external velocity in terms of its spatial 
Fourier transform. Since the fluid is incompressible, we will express the velocity in a 
general form that has no longitudinal part: 



where u is some general vector field. Next, we define the correlations of the velocity 
by the correlations of u, 

(ui(q, t)) = and 



where 5i m is the Kronecker delta, S() is the Dirac delta function and £ is the velocity 
correlation length. 




(1) 



(ui(q, ti)w m (p, h)) = 5i m 5(q + p)4>(£q)5{t 2 - h) 



(2) 
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5. The surface elements of the body are carried by the host fluid ||, i.e. each surface 
point moves according to 



r = v ext {r) +v^(f). 



(3) 



6. We assume that the external velocity is weak enough to cause only minor shape fluc- 
tuations of the body. 

We will be interested in the following in the Mean Squared Displacement (MSD) of the 
center. Since the body is deformable the definition of its center is not unique. We will 
choose a specific definition later. It turns out, however, that the value of the MSD at long 
times does not depend on the specific choice, because for long times the MSD (according to 
any definition) is much larger than the size of the body. 

Following the line of derivation of Schwartz and Edwards H0, equation ([5]) may be 
turned into a continuity equation for ■0, 



Consider a deformable body, carried by the host fluid in such a way that at any instant it is 
nearly spherical. Its state can thus be characterized by the position of its center, r*o(t), and 
a deformation function f(fl,t) that describes the shape by the equation 



where p is the distance of the surface from the center in the direction of the solid angle f2 and 
R is the radius of the body when not deformed. The deformation function / can be expanded 
in spherical harmonics, f(Q,t) = Y^oY,m=-i fl,m(t)Yi im (Q) . The center of the shape, r (t), 
is defined as that point around which fi m (t) = for m = —1, 0, 1. Parameterizing now the 
gauge field ift as ift(r, t) = + /(^> ~ 1; equation (f|) leads to a linear equation for each 



dip 
~dt 



+ (Vext + V<p) ■ V-0 = 0. 



(4) 



R 



+ f(n,t) -1 = 0, 



(5) 
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fl,mi 

9fi 1 

-jjjr + hfim + [p ■ (vext - r )]im = 0, (6) 

where p is a unit vector directed outwards from the center in the direction of Q, and 

[p ■ (v ext - ? )]im = Jdn{p- [v ext (f + R(l - f)p) - f ] Y* m (Q)} (7) 

The eigenvalues A/'s characterize the decay rate of a slightly deformed sphere into a sphere 
in the absence of the external velocity. Note that the term Xifi m in Eq. (|]) results from the 
shape induced velocity, in Eq. @) @]. Different physical systems are characterized by 
different sets of Aj's. For example, Schwartz and Edwards J7| calculate A; for a droplet with 
a constant surface tension and equal viscosities inside and outside the droplet, Gang at al. 
give A/ for a droplet with a surface tension under the assumption that the viscosity of 
its interior is much greater than that of the surrounding fluid, Milner and Safran || do the 
same for a surface controlled by bending energy and Dorries and Foltin add to it in-plane 
dissipation ||. In all of these cases, it is obvious that the decay must depend only on I 
because of the spherical symmetry. Therefore, our treatment is general and applicable to a 
large class of physical systems in which the decay rate depends only on I. 

Equation (||) implies that in order that /x, m stays zero for all times we must have as an 
equation determining the location of the center 

IP ■ {vext - r ))i m = ,m= -1,0,1. (8) 

For I ^ 1 it is clear that r*o can be dropped from the last term on the left hand side of Eq. (^j). 
Therefore f(£l, t) is linear in v ext (for long enough times the initial deformations have already 
decayed). Consequently we can always drop, for small enough v ext , f in the argument of 
v ext on the right hand side of Eq. (|7|). The result is decoupling of the deformation degrees 
of freedom from that of the center of the sphere. The equation for the motion of the center 
can thus be given, using linear combinations of Yi, m , in vector form as 

J dVt p (p • f ) = / dQ, p (p ■ v ext (f + Rp)) . (9) 



Direct integration of the left hand side, Fourier decomposition of v ex t(ro + Rp) on the right 



hand side and the use of the partial waves decomposition fllO| ,P 

CO I 



e 



-icf-(Rp) 



E E {-i) l ^Ji{qR)YL{Sl q )Y lm {V), (10) 



1=0 m=-l 

where fl q is the solid angle in the q direction, yields 

Fo = 3 1 dqe~^ (^UqR) + h(qR)A) v ext (q,t), (11) 
where jo an d J 2 are spherical Bessel functions and the matrix A(q) is given by 

A, = -lk+(k- q -f)- (12) 



It may seem that A on the right hand side of Eq. (|TT|) mixes directions. However, the 
bracketed term in Eq. ([E|) is just a projection operator on the transverse direction. The 
external velocity is incompressible and hence already transverse. Consequently, this term 
acts as a unity operator, 5ij and Eq. ( |TT| ) leads to 



?o = J dqe~^ ( Jo (qR) + j 2 (qR)) v ext (q, t). (13) 

Equation ( Jl3| ) is the explicit equation of motion for the center of the body. In the limit 
R — » the approximation, f = v ext (fo,t) is obtained. Note that this equation is general 
and describes the motion of the center for any given (small enough) external velocity field. 

Next, we calculate the MSD, ((Af r) function of the elapsed time, t. The starting 
point is the trivial relation 

Af (f)= f'rv(t')dt'. (14) 



Assuming the decomposition [12,13 



Wft>*i)W92>*2)e-^- w e-^< fa > ) = 

{v^MM) v extj (q 2 ,h)) (e-^We-^W) (15) 
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and using equations fll]), (@) and flT3|) the MSD is 

, . POO 

v (Af (t)) 2 ) = Stt J q 2 dq<P{q) (j (qR) + j 2 {qR)f t = 3Dt. (16) 



We expect the diffusion coefficient to depend on the ratio 7 = j in the following way. As 
7 increases the diffusion coefficient decreases. This is due to the fact that as 7 increases, 
different regions of the surface become less correlated and move in different directions. In the 
limit 7 — > 00, the movement of the center ceases and the MSD is always zero. In the limit 
7^0, the bracketed Bessel term in Eq. ( |L6j) can be replaced by unity. A close inspection of 
the derivation reveals that this limit produces the same MSD equation as the approximated 
equation r = v ext (ro) that is accurate in the limit of infinite correlation length, or point 
particles. 

We turn to evaluate the dependence of the diffusion coefficient on R and £ in the above 
limiting cases. Consider the correlation function <fr(q) = C (q£) a g(q£), where g is a cutoff 
function and g(0) = Cq > 0. The diffusion constant is given by 

87TC£ a f°° .,2+ Q /{„U. / \ , • / \12 



D = 3R^ Jo dU u 9(^u)\jo(u)+j 2 (u)} 2 . (17) 

In the limit R/£ —>■ 00 we distinguish between two cases: a < 1 and a > 1. Since the large 
u dependence of jo(u) +32(11) is proportional to cos{u)/u 2 we find that 



D oc <^ 



for a < 1 

R (18) 
for a > 1. 



In the opposite limit — > we find that regardless of a 



Doc^. (19) 



(Note here that we have written the power law dependence of 4>(q) as q a l; a but having other 
dimensional constants in the model may make C depend on £, so that eqs. flT%D and ( ^9] ) 
may be considered only as equations that yield the dependence of D on the radius R). 
It is interesting to consider the case where the fluctuations in the velocity field are due to 
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thermal agitation. General considerations show that the correlations of the external velocity 
must be, in this case, of the form 

{vextMM)v ex tM^)) = ( S *i ~ ^r) S(q + p)5(t 2 - h)-^, (20) 

where c is a dimensional constant (Strictly speaking the above is valid for qa > 1, where a 
is the inter-particle distance). We use these correlations with Eq. fll6|) and obtain 

DA (21) 

Dimensional analysis reveals that c must be proportional to KbT/t) (with a dimensionless 
proportionality constant). A detailed calculation yields a proportionality constant equal to 
(2ir)~ 3 . We do not represent here the detailed calculation, because it yields a result identical 



to that obtained in the past | I3] , |l5 | using two totally different approaches (that we wanted 
to check by a more direct calculation. Actually, Hadamard calculated the mobility that can 
be related to the diffusion constant by the Einstein relation). Therefore, for R larger than 
the inter-particle distance in the liquid, 

D = (22) 

Note that this result, for a liquid membrane that has liquid inside as well as outside, 
is different from the well known Stokes result for a hard sphere and from the result for a 
polymer subjected to thermal fluctuations |L6| (The difference is in the prefactors). 
Fig. Q depicts D for typical correlation functions while Fig. |2| depicts the dependence of 
the diffusion coefficient on R, £, for two noise realizations. The first is just a cutoff func- 
tion, <f)(q) = C exp(—q 2 ^ 2 ). The second corresponds to the temperature driven randomness, 
0(g) = C (q£)~ 2 exp(— q 2 ^ 2 )- We keep £ constant and vary R. There are two distinct regimes: 
For < 1 the diffusion coefficient is not sensitive to R, while for R/£ > 1 the slope of the 
graph turns towards -3 for the cutoff function and towards -1 for the temperature driven 
randomness. 
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Figure captions 



Fig. 1 

The dependence of the diffusion coefficient D on the ratio of the radius R of the sphere to 
the correlation length £, for typical correlation functions. The axis are presented in non- 
dimensional units, where -Do = C/£ 3 - 

Fig. 2 

The dependence of the diffusion coefficient on R. We keep £ constant and vary R. For 
negative ln(R/£) there is no R dependence, while for positive values D is proportional to 
R^ where \x tends to (a) -1 for thermal agitation and (b) -3 for a simple cutoff function. 
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FIG. 1. 
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i (a) M = C (qHr 2 e-^ 2 
I (b) </>{q) = C 
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FIG. 2. 



11 



